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(57) Abstract 

A digital radio communications svstem is presented in which the deleterious effects of fading, time dispersion and interference 
found in the channel are mitigated by using interference rejection and diversity combining. Digital signals representative of digital symbol 
sequences are transmitted through a radio channel in which the signals may be corrupted by fading, time dispersion and interference. The 
corrupted signals are received bv a receiver having one or more antennas. The received signals are processed to form a received signal 
sample stream r d (kT.) for each of the antennas. Estimated received signal samples for each of the antennas arc produced by combining 
hypothesized signal sequences and models of each channel. The received signal samples and the estimated received signal samples are 
processed in conjunction with an estimate of the impairment correlation properties to produce branch metrics to be used by a Viterbi. or 
other sequence estimation, algorithm. The branch metric processor may assume one of several configurations depending on whether or not 
the channel is time-varying. 
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METHOD OF AND APPARATUS FOR INTERFERENCE REJECTION COMBINING IN 
MULTI-ANTENNA DIGITAL CELLULAR COMMUNICATIONS SYSTEMS 

RELATED APPLICATIONS 

5 

This application is related to U.S. Patent Application Serial No. 

(attorney docket number E8699) entitled "Method of and Apparatus for Reducing the 
Complexity of a Diversity Combining and Sequence Estimation Receiver", by Johan 
Skold and Per-Olof Eriksson, and filed concurrently herewith. This related application 
10 is also assigned to the same assignee as the present application and is incorporated 
herein by reference in its entirety. 

FIELD OF THE INVENTION 

15 The present invention relates generally to the demodulation of digitally 

modulated radio signals received by a plurality of antennas and more specifically to the 
diversity combining of radio signals subjected to multipath fading, time-dispersion and 
interference. 

20 BACKGROUND OF THE INVENTION 

A common problem that occurs in the radio transmission of signals is that the 
signals are sometimes tost as a result of multipath fading and interference which may 
exist in a radio transmission channel. Throughout the following, the terms radio 

25 transmission channel, radio channel, and channel are used equivalent^ to refer to the 
same item. There are basically two multipath effects: flat fading and time dispersion. 
Rat fading arises from the interaction of the transmitted signal, or main ray. and 
reflections thereof, or echoes, which arrive at the receiver at approximately the same 
time. If there are a large number of reflections, flat fading gives rise to a Rayleigh 

30 distribution. Time dispersion occurs when the echoes are delayed with respect to the 
main ray. There may also exist in the radio environment signal sources which are not 
orthogonal to the desired signal. Non-orthogonal signals, or interference, often come 
from radios operating on the same frequency (i.e., co-channel interference) or from 
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radios operating on neighboring frequency bands (/.e., adjacent-channel Interference); 
non-orthogonal signal sources are referred to as interferes. 

One known method of reducing the influence of Raylelgh fading is to use a 
receiver having two or more mutually separated antennas, for instance as described in 
S Mobile Communications Oesign Fundamentals by William C. Y. Lee ( Howard W. 
Sams & Co., Indiana, USA. In section 3.5.1 of this book several examples are given 
describing how signals from two receiver amplifiers with separate antennas can be 
combined to counteract fading. These techniques are generally referred to as diversity 
combining. 

10 Time dispersion may be advantageously corrected by using an equalizer. In the 

case of digital signal modulation, a maximum likelihood sequence estimation (MLSE) 
equalizer such as described in Digital Communications, 2°* Ed., by John G. Proakis, 
Mc-Graw Hil! Book Company, New York, New York, USA, 1989 may be used. In section 
6.7 of this book, various methods are described for detecting signals corrupted by time 

15 dispersion, or inter-symbol interference (ISI), using MLSE equalization. 

The impact of other signal interference may be reduced by employing array 
processing techniques with multiple antennas. For example, adaptive beamforming can 
be used to "steer" a null in the antenna pattern in the direction of an interferes 

Recently, methods have been proposed that partially solve the problems of 

20 multipath fading and interference. In U.S. Patent 5,191,598 to Backstrom, et a/., for 
example, the problem of accurately detecting signals in the presence of fiat fading and 
time dispersion is overcome by using a Viterbi-algorithm having a transmission function 
estimated for each antenna. By reference thereto, U.S. Patent 5.191,598 is 
incorporated herein in its entirety by reference. 

25 A method of accurately detecting signals in the presence of flat fading and 

interference was presented in the IEEE Transactions on Vehicular Technology, Vol. 42, 
No. 4, Nov. 1993, J. H. Winters: "Signal Acquisition and Tracking with Adaptive Arrays 
in the Digital Mobile Radio System IS-54 with Flat Fading". 

In a practical radio communication system, there commonly co-exists flat fading, 

30 time dispersion, and interference. While the above-mentioned techniques address 
some of these problems, there exists a need to simultaneously and jointly combat flat 
fading, time dispersion, and interference. 
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Summary of the Invention 



In view of the foregoing background, it is therefore an object of the present 
invention to jointly combat and correct for the simultaneous existence of flat fading. 
5 time dispersion, and interference. 

A method is presented for generating and transmitting a signal representing a 
transmitted symbol sequence and receiving the signal on one or more antennas. The 
signal is processed to produce received signal samples for each of the antennas. 
Channel taps are estimated for each of the antennas. Impairment correlation properties 
10 are also estimated. Branch metrics are formed in a branch metric processor using the 
received signal samples and the channel tap and impairment correlation estimates. The 
branch metrics are employed in a sequence estimation algorithm to estimate the 
transmitted symbol sequence. 

In one embodiment, branch metrics are formed by generating hypothetical 
15 symbol sequences and filtering the hypothetical signal sequences with the channel tap 
estimates to produce hypothesized received signal samples for each antenna. The 
hypothesized received signal samples are subtracted from the received signal samples 
to produce hypothesized error signals which are processed with the estimate of 
impairment correlation properties to produce branch metrics. 
20 In another embodiment, an estimate of impairment correlation properties is 

formed by generating tentative detected symbol sequences which are filtered with the 
channel tap estimates to produce detected signal samples lor each antenna. The 
detected signal samples are subtracted from the received signal samples to produce 
detected error signals which are processed with the estimate of impairment correlation 
25 properties to produce an update of the estimate of the impairment correlation 
properties. 

In still another embodiment, branch metrics are performed by generating 
hypothetical symbol sequences and calculating pre-computed values for all of the 
hypothetical symbol sequences using the channel tap estimates and the estimates of 
30 impairment correlation properties. The pre-computed values are processed with the 
received signal samples to produce branch metrics. 

In yet a further embodiment branch metrics are formed by generating 
hypothetical symbol sequences and combining the received signal samples with the 
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channel tap estimates and the estimate of impairment correlation properties to produce 
metric multipliers which are processed with the hypothetical symbol sequences to 
produce branch metrics. 

An apparatus is also presented comprising a digital transmitter for generating 
5 and transmitting a signal representing a transmitted symbol sequence, a receiver for 
receiving said signal on one or more antennas, an analog to digital converter for 
converting the signal into received signal samples for each of the antennas, a channel 
tap estimator for estimating channel taps for each of the antennas to produce channel 
tap estimates, an impairment correlation estimator for estimating impairment 

10 correlation properties to produce an estimate of impairment correlation properties. The 
apparatus also comprises a branch metric processor coupled to the channel tap 
estimators, the impairment correlation estimator and the received signal samples for 
forming branch metrics using the received signal samples, the channel tap estimates, 
and the estimate of impairment correlation properties. The output of the branch metric 

15 processor is coupled to a sequence estimation processor which produces an estimate 
of said transmitted symbol sequence. 

In another embodiment the branch metric processor comprises a symbol 
sequence generator for generating hypothetical symbol sequences coupled to a digital 
filter for filtering the hypothetical signal sequences with the channel tap estimates to 

20 produce hypothesized received signal samples for each antenna. The filters are 
coupled to means for subtracting said hypothesized received signal samples from the 
received signal samples to produce hypothesized error signals which are coupled to 
means for processing said hypothesized error signals with said estimate of impairment 
correlation properties to produce branch metrics. 

25 Yet another embodiment is presented in which the impairment correlation 

estimator comprises a symbol sequence generator for generating tentative detected 
symbol sequences coupled to a digital filter for filtering the tentative detected symbol 
sequences with the channel tap estimates to produce detected signal samples for each 
antenna. The filter outputs are coupled to means for subtracting said detected signal 

30 samples from the received signal samples to produce detected error signals which are 
coupled to means for processing said detected error signals with said estimate of 
impairment correlation properties to produce an update of the estimate of impairment 
correlation properties. 
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These and other features and advantages of the present invention will be 
readily apparent to one of ordinary skill in the art from the following written description 
when read In conjunction with the drawings In which like reference numerals refer to 
like elements. 

5 Brief Description of the Drawings 

Figure 1 is a schematic illustration of a digital radio communication system; 

Figure 2 is a schematic illustration of a receiver processor and transmission function 
10 according to the present invention; 

Figure 3 is a schematic illustration of a transmission function; 

Figure 4 is a schematic illustration of a branch metric processor according to the 
15 present invention; 

Figure 5 is a schematic illustration of an adaptive estimator of the impairment 
correlation properties; 

20 Figure 6 is a schematic illustration of yet another embodiment of a branch metric 
processor according to the present invention; 

Figure 7 is a schematic illustration of another embodiment of a receiver processor and 
a transmission function according to the present invention; and 

25 

Figure 8 is a schematic illustration of yet another embodiment of a receiver processor 
and a transmission function according to the present invention. 

DESCRIPTION OF THE INVENTION 

30 

In the following description, for purposes of explanation and not limitation, 
specific details are set forth, such as particular circuits, circuit components, techniques, 



'.SOOCIO: <WO 96047 38A1 J_> 



WO 96/04738 



PCT/US95/09484 



etc- In order to provide a thorough understanding of the invention. However It will be 
apparent to one of ordinary skill In the art that the present invention may be practiced in 
other embodiments that depart from these specific details. In other Instances, detailed 
descriptions of well-known methods, devices, and circuits are omitted so as not to 

5 obscure the description of the present invention with unnecessary detail. 

A radio transmitter and receiver system for a radio communication system is 
illustrated schematically in Figure 1. The radio communication system may operate 
using frequency division multiple access (FDMA). time division multiple access (TDMA), 
or code division multiple access (CDMA), or some combination thereof. A transmitter 

10 has a digital symbol generator 102 which receives an information carrying signal 101 
and generates a corresponding digital symbol sequence, S. The symbols S are 
subjected to digital to analog (D/A) conversion, modulation, pulse shape filtering, 
amplification, and are transmitted as analog signal Y by digital transmitter 103 
according to known techniques. 

15 In addition to thermal noise, there may also exist an interferer 108 transmitting 

signal X which may be non-orthogonal to signal Y. Signals Y and X travel through 
separate radio channels and are intercepted by antennas 104 which are D in number. 

Radio units 105 amplify, downconvert, and filter the received signals according 
to known methods to produce analog signals. Each analog signal is coupled to an 

20 analog-to-digital (A/D) converter 106, which converts the analog signal into a received 
signal sample stream r d {kT t ), where T t is the sample period, the reference numeral k is 
an integer counter, and the subscript d indicates that the signal arrives from the d m 
antenna 1 £ d S D. The sampling period T s may be less than the symbol period T. The 
received signal sample streams are collected in processor 107, which processes these 

25 streams to produce an estimate of the transmitted digital symbol stream, § . In later 
descriptions, transmission function 109 is used to refer to the signal path through digital 
transmitter 103, the radio transmission channel (not shown in Figure 1), antennas 104, 
radio units 105 and A/Ds 106 collectively. 

The processing unit 107 is illustrated in greater detail in Figure 2 where, for 

30 simplicity, the number D of antennas is restricted to three: designated a, b, and c. 
Processing unit 107 may be, for example, a Digital Signal Processor (DSP) such as a 
TMS320C50 manufactured by Texas Instruments. The function of processing unit 107 
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is to produce an estimate of the transmitted digital symbol stream S which accurately 
corresponds to the symbol sequence S which was originally transmitted. 

Transmission function 109 produces the received signal sample streams r t (kT,) ( 
r b (kT,) ( and r e (kT f ) which are sent to processing unit 107 where they are processed in 
5 accordance with the present invention. The received signal sample streams r t (kT $ ), 
rb(kT,) t and r c (kT,) are coupled to a signal pre-processor, or sync, block 206 where the 
received signal sample streams are correlated with known timing/synchronization 
sequences as described, for example, by Giovanna, etat. "Fast Adaptive Equalizers for 
Narrow-Band TDMA Mobile Radio", IEEE Transactions on Vehicular Technology, Vol. 
10 40, No. 2, May 1991 , pp. 392-404. For the case of symbol-spaced demodulation, if the 
sample period T t is less than symbol period, T, the signal pre-processor 206 performs a 
decimation of the received signal sample streams r.(kT t ), r b {kT t ), and r c (kT $ ) to produce 
one sample per symbol: designated r.(n), r b (n), and r c (n) respectively. For the case of 
fractionally-spaced demodulation, more than one sample per symbol is generated. 
15 Estimating circuits 202a, 202c, and 202c produce channel tap estimates cCc), 

Ct>(T), and Cc(x) which are used to model the radio transmission channel associated with 
each particular antenna. Initial channel tap estimates can be obtained from sync 
correlation values or least-squares estimation according to known techniques. If the 
channel must be tracked, it is typical to use received data and tentative symbol 
20 estimate values generated in the sequence estimation processor 204. Channel 
tracking is known to those skilled in the art as discussed, for example, in Digital 
Communications 2nd Ed. by Proakis as previously mentioned, and by A. P. Clark and 
S. Hariharan, 'Adaptive Channel Estimates for an HF Radio Link - , IEEE Trans, on 
Communications, vol. 37, pp. 918-926, Sept. 1989. The channel tap estimates c(t). 
25 Ct(-c), and Cc(t) are coupled to the input of the branch metric processor 203. 

Also coupled to the branch metric processor 203 is an estimate of the 
impairment correlation properties obtained from impairment correlation estimator 207. 
The estimate of the impairment correlation properties comprises information regarding 
the impairment correlation properties between the antennas 104, or between relative 
30 sample phases as will be discussed in greater detail hereinafter. The impairment 
correlation estimator uses impairment process estimates to update and possibly track 
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the estimate of the impairment correlation properties which is discussed in further detail 
in the ensuing text and figures. 

Branch metric processor 203 uses received signal samples r,(n), r 6 (n), and r e (n) ( 
channel tap estimates c t (t) ( Ct>(*c) ( and Cc(t) f and the estimate of the Impairment 
correlation properties to form branch metric Mh(n). This branch metric is used, for 
example, In a sequence estimation processor 204 to develop tentative and final 
estimates of the transmitted symbols. 

It will be appreciated by one of ordinary skill in the art how the present invention 
operates when fractionally-spaced (T/M) equalization is employed. In this case, each 
signal preprocessor 206 produces M samples per symbol period, which correpond to 
M different sampling phases. This gives rise to M received signal sample streams per 
antenna. As a result, branch metric processor 203 receives D times M sample streams. 
As before, the channel tap estimators 202a - 202c produce channel tap estimates for 
antennas, a • c respectively. However, these channel tap estimates are more in 
number, as estimates for each sampling phase must be provided. AJso, the impairment 
correlation estimator 207 estimates impairment correlation properties over antennas 
and sampling phases. This is done by treating the D antennas and M sampling phases 
as D times M antennas. As a result, the branch metric processor operates in a manner 
equivalent to having D times M antennas and one sample per symbol. It is apparent 
that the present invention can be applied to the case of one receive antenna and 
fractionally-spaced equalization. 

The transmission function 109 is illustrated in greater detail in Figure 3 where, 
for simplicity, the number of interferers is restricted to one. It is obvious to one skilled in 
the art that the present invention may also be used for the case where there are two or 
more interferers. The transmit function 109 begins with the signal path for the symbol 
sequence S through digital transmitter 103 which transmits analog signal Y. The analog 
signal Y propagates through a separate radio transmission channel to each of the three 
receiver antennas: radio channel 301a to receiver antenna 104a, radio channel 301b to 
receiver antenna 104b, and radio channel 301c to receiver antenna 104c. Similarty, 
interference signal X also propagates through three other separate radio channels 
302a-302c to receiver antennas 104a-104c respectively. Radio channels 301a-301c 
and 302a-302c may introduce fading and time dispersion. Omnipresent thermal noise 
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processes n, • n e are also received by receiver antennas 104a-104c respectively. Each 
antenna 104a«104c Is coupled to a radio unit 105a-105c respectively which amplifies, 
downconverts, and filters the received signals according to known methods to produce 
an analog signal. Each analog signal is coupled to an analog-to-dlgital (A/D) converter 
5 I06a-106c which converts the analog signals into received signal sample streams 
r»(kT,). r b (kT,), and r e (kT t ). One method for conversion from analog to digital is to use 
log-polar signal processing, as described in U. S. Patent 5,048.059 to Dent 
incorporated by reference herein in its entirety. For subsequent processing, a 
conversion from log-polar to rectangular samples is made, so that, for example, I and 
10 Q samples, sometimes referred to as complex samples, are used. By using log-polar 
signal processing initially, a limiting receiver which provides signal strength and phase 
samples can be used, and adaptive gain control can be made simple. 

In an MLSE equalizer, all possible transmitted symbol sequences S are 
considered. In one implementation, hypothesized symbol values s*(n) are filtered by 
15 channel tap estimates c.(t), c(t). and Cc(t) to produce hypothesized received samples 
r.h(n), r b . h (n), and r e .„(n) for each antenna. The differences between the hypothesized 
r M (n)-r e . h (n) and the actual r.(n)-r e (n) received signal sample streams, referred to as the 
hypothesis errors, give an indication of how good a particular hypothesis is. The 
squared magnitude of the hypothesis error is used as a metric to evaluate a particular 
20 hypothesis. The metric is accumulated for different hypotheses for use in determining 
which hypotheses are better using the sequence estimation algorithm. This process 
may be efficiently realized using the Viterbi algorithm, which is a known form of 
dynamic programming. A description of the Viterbi alrogithm may be found in Forney, 
G.. The Viterbi Algorithm', Proc. of the IEEE, vol. 61 , pp. 268-278. March. 1 973. As is 
25 obvious to one of ordinary skill in the art, other sequence estimation algorithms, such 
as the M-algorithm. may also be used. 

In an MLSE equalizer, there are states associated with different symbol 
sequence hypotheses s„(n). At a given iteration, there are previous states: each 
associated with an accumulated metric. Each pairing of a previous state with a current 
30 state results in a branch metric. M h (n). The candidate metric for a current state is then 
the sum of the branch metric M h {n) and the previously accumulated metric. For each 
current state, the previous state which gives the smallest candidate metric is selected 
as the predecessor state, and the smallest candidate metric becomes the accumulated 
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metric for the current state. For metric combining, as described in aforementioned U.S. 
Patent 5,191,598, the branch metric can be expressed as: 

M h (n) = [r(n) - Cs,(n)f D[r(n) - Ca h (n)] 
where: 

r(n) = [r g (n)r b (n)r c (n)] 



C = 



c a (0) 

c b (0) 
c c (0) 



c 4 (N t -1) 
c to (N t -1) 
c c (N t -1) 



s h (n) = [s h (n)s h (n-1)...] T 



10 



K, 0 0 

0 K b 0 
0 0 K c 



The channel tap estimates for each signal received on antenna 104a-104c are 
designated by c(t), (^(x), Cg(-c) respectively where t is the delay (i.e. x - 0 is the main 
ray, t = 1 is the first echo, etc.). N, is the number of channel taps estimated per antenna 
and K«, K*. Kc are weighting coefficients for antennas 104a- 104c respectively. 

The present invention takes advantage of the fact that, from a diversity and 
equalization point-of-view t the impairment (interference + noise) on multiple receive 
antennas 104 is often correlated. By expanding diversity combining techniques to 
exploit this correlation, significant gains are realized. For optimal performance, a 
whitening, or decorrelation, process may be applied and the optimum branch metric 
may include the inverse of the impairment correlation matrix. The optimum branch 
metric M h (n) according to the present invention is: 
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M h (n) = [r(n) - C(n)t h (n)] H A(n)[r(n) - C(n)a h (n>] = # h M (n) A(n)« h (n) 
where: 

A(n) = R a (n)-', or a related quantity; R n (n) = E(z(n)z H (n)) 
2(n) = [z, (n)z 6 (n)z e (n)] T ; e ft (n) = r(n) - C(n) s„ (n) 

The time varying nature of the channel and the impairment correlation are denoted with 
time index n. The R B (n) matrix is referred to as the impairment correlation matrix at 
discrete time, n. The A(n) matrix (i.e., the A-matrix) is the inverse of the R n (n) matrix, or 
a related quantity such as the adjoint or pseudo-inverse. As is obvious to one of 
ordinary skill in the art, R a (n) and A(n) are specific examples of impairment correlation 
properties of which other forms are known. Throughout the following, the term A-matrix 
is used genericaily to refer to any estimate of the impairment correlation properties. 

The impairments on antennas 104a-104c at time n are designated by z.(n), 
z,,(n), and Zc(n) respectively. For a given hypothesis, e h (n) is an estimate of the 
impairment process. As shown above, the A-matrix. A(n). is the inverse of the 
impairment correlation matrix R a (n). For the case of uncorrelated impairment (i.e., no 
interferer) the A-matrix reduces to diagonal matrix D. When the signal is known or 
detected correctly, the impairment is given by: 

z(n) = r(n)-C(n)s d .,(n) 

15 where: 



10 



•d..(n) = [s a .,(n)s d .,(n-1)...] T 
Note that s„«(n) is the known or detected symbol sequence at time n. 

Determination of the A-matrix for use in the present invention can be performed 
in a number of ways depending upon the specific application and the required 
performance. The simplest approach is to use a fixed set of values for the A-matrix, 
20 stored in memory, that are never updated. These values depend primarily on the 
configuration of the receive antennas and on the carrier frequencies being employed. 
An alternative approach is to determine the A-matrix from synchronization information 
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and to keep the A-matrix values constant between synchronization, or other, known 
fields. At each new occurrence of the synchronization field, the A-matrix can be 
recomputed, with or without use of the previous A-matrix values. Another alternative 
approach is to use synchronization fields to initialize, or improve, the A-matrix values 

5 and then to use decisions made on the data field symbols to track the A-matrix values. 

Also, consideration is given for the method used to track the A-matrix values. 
Since the A-matrix comprises information regarding the impairment correlation 
properties standard estimation methods for estimating correlation or inverse correlation 
matrices can be applied. Using either known or detected symbol values, impairment 

10 values can be obtained by taking the differences between the received signal sample 
streams r t (n)-r c (n) and the hypothesized received signal sample streams r t . h (n)-r Cth (n). 
At time n, this gives a vector of impairment values, denoted z(n); one value for each 
antenna, A straightforward way of forming the A-matrix is given by: 

fl(n) = XR(n-1) + Kz(n)z H (n) 

A(n) = FT 1 (n) 



15 K is a scaling constant, typically 1 or . Because R(n) is a Hermitian matrix, only 

a portion of the matrix elements need be computed. 

Such a straightforward approach is fairly high in complexity. One way to reduce 
complexity is to apply the matrix inversion lemma and update the A-matrix directly as: 



where: 



A < n - 1) - 1 7 — tVtt I p<") p h < n > 

1 X + z(n) M p(n) 1 



p(n) = A(n-1)2(n) 

20 Because the A-matrix is Hermitian, it is only necessary to compute those elements on 
the diagonal and either those elements above or below the diagonal. 

These techniques for estimating and tracking the A-matrix are given only for 
purposes of illustration. In general, the A-matrix can be expressed and estimated in a 
variety of ways, as is obvious to one of ordinary skill in the art. See. for example, the 

25 book by S. Haykin, Adaptive Filter Theory, Second Edition, Prentice-Hall, Englewood 
Cliffs, N.J., 1991- The present invention may also be applied to the blind equalization 
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problem In which known synchronization sequences are absent. In this case, the A- 
matrix is estimated in a manner similar to how the channel is estimated. 

For the purposes of illustration, the present invention will now be described In 
greater detail as exemplified in four different embodiments. 



Description of Embodiment having a Symbol-Spaced Equalizer with Two Antennas 
and Adaptive Channel Estimates 

10 

In a first embodiment, a processor 107 having a symbol-spaced (i.e., T-spaced) 
equalizer is presented where the channel must be tracked over the data field, or burst. 
This embodiment is applicable to the digital cellular system defined by the IS-54B 
specification which has relatively long TDMA data bursts, or time slots, (6.67 

15 milliseconds) with respect to time. For this embodiment, the branch metric processor 
203 is illustrated in greater detail in Figure 4 where, for simplicity, the number of 
antennas is further restricted to two: designated a and b. This particular embodiment 
has usefulness in that the use of two receive antennas is common in many cellular 
systems which already employ some form of diversity combining. As before, it is 

20 obvious to one of ordinary skill in the art that this embodiment may also be employed in 
the case where there are three or more antennas. 

The impairment correlation matrix R u and the inverse impairment correlation 
matrix A are defined as follows: 



' Pt* 


-p*" 


= w 


' Pm, 


-P*' 




Pm . 









25 The variable p„ denotes the impairment power received on antenna a; the variable p*, 
denotes the impairment power received on antenna b. The off-diagonal matrix 
elements are the cross correlation values: denotes the correlation of the impairment 
received on antenna a with the conjugate of that received on antenna b. 
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The branch metric then becomes: 
M h (n) = e h (n) H A(n^^ 

where: 

PmP» -|Pib| 

The calculation of this branch metric is schematically illustrated in Figure 4. 

A symbol sequence generator 410 generates hypothesized symbol sequences 
5 s^n). These sequences are filtered in filters 400 using channel tap estimates c a (t) and 
Ct,(x) for antennas a and b to produce hypothesized received signal samples r ih (n) and 
ro.h(n), respectively. The hypothesized received signal samples, r M (n), are subtracted 
from the actual received signal samples from antenna a ( r,(n), in summing junction 401 
to produce error signals e^(n). Similarly hypothesized received signal samples, ^(n), 

10 are subtracted from the actual received signal samples from antenna b» r b (n), in 
summing junction 402 to produce error signals e b .„(n). Blocks 403 form the squared 
magnitudes of the error signals e^fn) and e b(h (n). The squared magnitude for error 
signals e.. h (n) is multiplied at junction 406 by multiplier m Mt the result being coupled to 
summing junction 408. The squared magnitude for error signals e b ^(n) are multiplied at 

15 junction 407 by multiplier m*,, the result being coupled to summing junction 408. 
Finally, multiplier 404 forms the product of e a . h (n) and e\ h (n) t the product of which is 
subsequently multiplied by multiplier m* in multiplier 405, forming the real part only. 
The result is subtracted in summing junction 408, the output of which is the branch 
metric Mh(n). The multipliers m M , m^, and m^ are related to the impairment correlation 

20 matrix by: 

m M = wptt 

m at> = Swp^ 

As is obvious to one skilled in the art the w term is common to the branch metric 
calculation and may be applied in a different manner or even omitted when the 
denominator to w approaches zero. 
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At time n, the A matrix elements are updated as: 
p tt (n + 1)=Xp tt (n) + |e.(n]( J K 

p* (n + 1) = Xp^ (n) + e 4 (n)e b ' (n)K 

P tt> (n + 1)=Xp be (n) + |e b (n]j 2 K 

K is a scaling factor which, if equal to unity, is dropped from the calculation to reduce 
the number of computations. K may be derived from X which is the so-called "forgetting 
5 factor*. 

A schematic illustration of the impairment correlation matrix update Is shown in 
Figure 5. Tentative detected symbol values s^n) from sequence estimation processor 
204 are filtered in filters 500 using channel tap estimates c,(t) and Cb(T) from channel 
tap estimators 202 for antennas a and b to produce expected received samples r^fn) 

10 and r b .det(n) ( respectively. Impairment signal z 4 (n) is produced by subtracting, in 
summing junction 501 r tid «(n) from the actual received signal samples from antenna a, 
r t (n). Similarly, impairment signal z*(n) is produced by subtracting, in summing junction 
502 r bi «*,(n) from the actual received signal samples on antenna b, r b (n). If the tentative 
detected symbol values are correct and the channel tap estimates are accurate, then 

15 error signals z«(n) and Zt>(n) represent the impairment received on antennas a and b 
respectively. Impairment signals z«(n) and z^n) are scaled by the root of the scaling 
factor K in multipliers 503 and 505 respectively to produce scaled impairment signals 
which are coupled to blocks 506 and 507, respectively. 

The impairment power received on antenna a, p,i(n), is multiplied in multiplier 

20 511 by the forgetting factor, X, and summed in junction 510 with the squared magnitude 
of the scaled impairment signal from block 506 to produce the updated impairment 
power p M (n+1). The value of p M (n+1) is then used to overwrite the memory location 
515 of the previous impairment power p M (n). Similarly, the previous impairment power 
received on antenna b, Pw>(n), is multiplied in multiplier 513 by the forgetting factor, X, 

25 and summed in junction 512 with the squared magnitude of the scaled error signal from 
block 507 to produce the updated impairment power Pbt>(n+1) which is used to overwrite 
the memory location 514 of the previous impairment power p«>(n). To produce the 
updated impairment cross correlation, the scaled error signal from multiplier 503 is 
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multiplied with the conjugate of the scaled error signal from multiplier 505 In junction 
504. Also, the previous cross-correlation p^n), stored in memory 51 s, Is scaled by the 
forgetting factor in multiplier 509. The output of junction 504 Is summed In junction 508 
with the output of multiplier 509 to yield the updated cross correlation p*(n+1). As 
5 before, the updated value p^n+l) is used to overwrite the memory location 516 of the 
previous value Pi*>(n). 

There is typically a delay in updating the channel tap estimates which allows the 
tentative detected symbols to become reliable. In U.S. patent 5,164,961 by 
Gudmundson, et a/., incorporated herein by reference in its entirety, this delay is 
10 avoided by using multiple channel models: one for each state in the sequence 
estimation processor 204. With the present invention, there is also a delay in updating 
the A-matrix quantities. It is obvious to one of ordinary skill in the art that this delay can 
be avoided by using multiple A-matrices; one for each state in the sequence estimation 
processor 204. 

15 

Description of Embodiment having a Symbol-Spaced Equalizer with Two Antennas 
and Fixed Channel Estimates 

In a second embodiment a receiver having a symbol-spaced equalizer is 
20 presented where the channel can be considered static over a TDMA burst. This 
embodiment is applicable to the digital cellular system defined by the pan-European 
GSM specification which has relatively short (577 microseconds) TDMA data fields with 
respect to time. In this case, the impairment correlation matrix R n , and hence the A- 
matrix, do not vary over the burst and can be estimated from the sync word embedded 
25 in the GSM frame. In this instance, the branch metric can be expressed as: 
M h (n) = [m - Cs h (n)f A(r(n) - Cs h (n)) 

M h (n) = r H (n) Ar<n) - 2Re{s h (n) H C H Ar<n)} + s h (n) H C H ACs h (n) 

Referring now to Figure 6, processing time is saved by pre-computing and storing in 
memory 601 the following values for all possible hypothetical symbol sequences, a*: 
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f h = 2A H Ct h 
g h =:S h H C H ACs h 

Using these pre-computed values, the branch metric may be simplified as: 

M h (n) = -Re{f h H r(n)} + g h 

wherec 

r(n) = [r a (n)r b (n)] T 

For each hypothesized symbol sequence 8 h (n), the received signal samples r a (n) and 
5 r b (n) are multiplied in blocks 602 and 603, respectively, by the corresponding f h value, . 
which was precomputed as described above and retrieved from memory 601; only the 
real parts of the products are formed. The h-index refers to the hypothesis index which 
relates a particular f h value to a particular hypothesized symbol sequence a*(n). The 
output of multipliers 602 and 603 are coupled to summing junction 604 where they are 
10 summed and the results passed to junction 605 where the output of junction 604 is 
subtracted from the corresponding g h value: also precomputed as described above and 
retrieved from memory 601. The result is the branch metric M h (n). 



Description of Embodiment having a Fractionally-Spaced Equalizer with one or 
15 more Antennas and Adaptive Channel Estimates with Partial Optimization 



In a further embodiment a receiver having one or more antennas coupled to a 
processor 107 comprising a fractionally-spaced (TfM) equalizer with adaptive channel 
tracking is presented. This particular embodiment is applicable to the digital cellular 

20 system defined by the IS-54B specification which has relatively long (6.67 milliseconds) 
time slots, or bursts, such that the channel tap estimates need to be updated during 
the duration of the burst. For this embodiment, receiver processor 107 is illustrated in 
greater detail in Figure 7 where, for simplicity, only a single interferer X and three 
antennas designated a, b, c are considered and M=2. It is obvious to one of ordinary 

25 skill in the art that the present invention may be practiced with a plurality of interferes 
and a plurality of antennas and a fractional spacing other than M=2. The transmission 
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function 109 shown in Figure 7 is, for example, Identical to that described by Figure 3 
which produces three received signal sample streams: r t (kT f ), r b (kT a ) and r c (kT,). The 
received signal sample streams: r t (kT t )-r c (kT t ) are each coupled to a signal 
preprocessor 707 which performs correlations with a known synchronization pattern 

5 similar to that described for block 206 of Figure 2. Timing is determined from the 
correlation values using some optimization criterion, such as maximum energy in the 
channel taps. In this particular embodiment, since M has been chosen equal to 2, each 
signal pre-processor 707 produces, from each incoming received signal sample stream, 
two samples per symbol. For example, as shown in Figure 7, r a , 0 (n), r tt1 (n) are 

10 generated from received signal sample stream r^kT,). Similarly r bi0 (n), r bt1 (n) and rco(n) ( 
rct(n) are produced from received signal sample streams r b (kT t ) and r c (kT t ) 
respectively. If, for example, M is chosen to be equal to 4 then 4 samples per symbol 
would be produced. These signals, r^n), r^n), r bi0 (n), r bt1 (n), r^n), and ^(n), are 
coupled to the branch metric preprocessor 701 where they are treated as though each 

15 signal, r tt0 (n), r,.,(n), r bi0 (n), r bt ,(n) r c0 (n), and r CJ (n), arrived from a separate antenna. 

In practice, one can expand the expression for the branch metric and collect 
computations to provide a lower complexity implementation. Two techniques can be 
used to reduce complexity: a) expand the branch metric expression and collect like 
terms, and b) reorder at which iterations certain terms are computed. Using the first 

20 technique, metric multipliers, which can be pre-computed for all hypotheses, are 
expressed as: 

e(j,n) = c H (j,n)A<n)r(n) 

fan) = c H an)A(n)c(j,n) 

gG,M = e H G»n)A<n)c<k t n) wherek>j 
where cQ.n) is the jth column of C(n). J is the number of channel taps. In other words, 

the received signal vector r(n) is modeled as: 

j 

25 r(n) = ^s h (n-j)c(j,n) 

ro 

The values eQ.n), f(j,n), and g(j.k,n) are computed by the branch metric preprocessor 
701 and stored in memory (not shown). The indices j and k are ray indices. The values 



8NSDOCID: <WO 9604738A1J_> 



WO 96/04738 



PCT/US95/09484 



19 

e(],n), f(j,n), and g(J.M) are referred to as metric multipliers since they can be used in 
multiplications to form the equalizer metrics in branch metric processor 704. The values 
e(].n), f(j.n), and g(],k,n) are coupled to branch metric processor 704 which computes 
the branch metric M h (n) according to: 

M h (n) = -X 2Re { e(i '^ 
5 W ° h ° * 

+2 L £Re{9(^ k ' n ) s ^i n -l>^ n " k )} 

In systems where all symbols have the same amplitude, such as binary phase 
shift keying (BPSK) and quaternary phase shift keying (QPSK) systems, the term 
| Sh (n-T)| 2 is a constant, independent of hypothesis. Thus, this term can be dropped, 
obviating the need to compute and use the f(j,n) metric multipliers. 
10 While the embodiment was illustrated for the case of three antennas, it is 

obvious to one of ordinary skill in the art that the present invention may be practiced 
with a single antenna and fractionally-spaced equalization. 

Description of Embodiment having a Fractionally-Spaced Equalizer with one or 
15 more Antennas and Adaptive Channel Estimates with Full Optimization 

A further reduction in complexity can be achieved by employing the receiver 
processor 107 shown in Figure 8. The receiver processor 107 of Figure 8 is essentially 
identical to that of Figure 7 with the difference being the branch metric preprocessor 
20 801 and the branch metric processor 804. As shown in Figure 8, branch metric 
preprocessor 801 pre-computes and stores in memory (not shown) the following 
alternate metric multipliers: 

Z(n) = X c " + j) ^ n + rtr(n + & 

i 

S(k, n) = £c H (j,n + j)A(n + j)c(j + k,n + j) 

In practice, it may be advantageous to approximate A{n+j) with A(n) and c(j,n+j), 
25 c(j+k,n+j) with cQ.n), c(j+k,n) respectively. Alternate metric multipliers Z(n) and S(k,n) 
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are coupled to branch metric processor 804 which computes the branch metric 
according to: 

M h (n) = Re s h ' (n) 2Z(n) + s h (n)S(O.n) + J 2S(k, n)s h (n - k) 
L L hii 

While the embodiment was illustrated for the case of three antennas, it is 
5 obvious to one of ordinary skill in the art that the present invention may be practiced 
with a single antenna and fractionally-spaced equalization. While the present invention 
has further been described with respect to a particular digital cellular communications 
system, those skilled in the art will recognize that the present invention is also 
applicable to other communications systems and that therefore the present invention is 

10 not limited to the specific embodiments described and illustrated herein. Different 
embodiments and adaptations besides those shown and described as well as many 
variations, modifications and equivalent arrangements will now be reasonably 
suggested by the foregoing specification and drawings without departing from the 
substance or scope of the invention. While the present invention has been described 

15 herein in detail in relation to its preferred embodiments, it is to be understood that this 
disclosure is only illustrative and exemplary of the present invention and is merely for 
the purposes of providing a full and enabling disclosure of the invention. Accordingly, rt 
is intended that the invention be limited only by the spirit and scope of the claims 
appended hereto. 
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What is claimed is: 



1 1. A method for reducing the effects of signal fading, time dispersion, and Interference 

2 in a radio communications system, comprising the steps of: 

3 (a) generating and transmitting a signal representing a transmitted symbol 

4 sequence; 

5 (b) receiving said signal on at least two separate antennas; 

6 (c) processing said signal to produce received signal samples for each of said 

7 antennas; 

8 (d) estimating channel taps for each of said antennas to produce channel tap 

9 estimates; 

10 (e) estimating impairment correlation properties to produce an estimate of 

11 impairment correlation properties; 

12 (f) forming branch metrics in a branch metric processor using said received 

13 signal samples, said channel tap estimates, and said estimate of impairment correlation 

14 properties; and 

15 (g) employing said branch metrics in a sequence estimation algorithm to 

16 estimate said transmitted symbol sequence. 

1 2. A method according to claim 1, in which processing step (c) employs log-polar signal 

2 processing. 

1 3. A method according to claim 1, in which said estimate of impairment correlation 

2 properties is represented as the inverse of an impairment correlation matrix. 
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1 4. A method according to claim 1 , in which said forming step (f) comprises: 

2 (a) generating hypothetical symbol sequences; 

3 (b) filtering said hypothetical signal sequences with said channel tap estimates 

4 to produce hypothesized received signal samples for each antenna; 

5 (c) subtracting said hypothesized received signal samples from the received 

6 signal samples to produce hypothesized error signals; and 

7 (d) processing said hypothesized error signals with said estimate of impairment 

8 correlation properties to produce branch metrics. 

1 5. A method according to claim 1 , in which said estimating step (e) comprises: 

2 (a) generating tentative detected symbol sequences; 

3 (b) filtering said tentative detected symbol sequences with said channel tap 

4 estimates to produce detected signal samples for each antenna; 

5 (c) subtracting said detected signal samples from the received signal samples to 

6 produce detected error signals; and 

7 (d) processing said detected error signals with said estimate of impairment 

8 correlation properties to produce an updated estimate of said impairment correlation 

9 properties. 

1 6. The method of claim 1 , in which said forming step (0 comprises: 

2 a) generating hypothetical symbol sequences; 

3 b) calculating pre-computed values for all of said hypothetical symbol 

4 sequences using said channel tap estimates and said estimate of impairment 

5 correlation properties; and 

6 c) processing said received signal samples with said pre-computed values to 

7 produce branch metrics. 
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1 7, A method according to claim 1 , In which said forming step (f) comprises: 

2 a) generating hypothetical symbol sequences; 

3 b) combining the received signal samples with channel tap estimates and said 

4 estimate of Impairment correlation properties to produce metric multipliers e(j,n); 

5 c) combining said channel tap estimates with one another and said estimate of 

6 impairment correlation properties to produce metric multipliers f(j,n) and g(j,k t n); and 

7 d) combining said metric premultipliers with said hypothetical symbol sequences 

8 to produce branch metrics. 



1 8. A method according to claim 1 , in which said forming step (f) comprises: ; 

2 a) generating hypothetical symbol sequences; 

3 b) combining the received signal samples with channel tap estimates and said 

4 estimate of impairment correlation properties to produce metric multipliers e(j,n); 

5 c) combining said channel tap estimates with one another and said estimate ,df 

6 impairment correlation properties to produce metric multipliers g(j,k,n); and 

7 d) combining said metric premultipliers with said hypothetical symbol sequences 

8 to produce branch metrics. 

9 



1 9. The method of claim 1, in which said employing step (g) further comprises the use of 

2 the Viterbi algorithm. 

1 10. A method according to claim 1, in which said forming step (f) comprises: 

2 a) generating hypothetical symbol sequences; 

3 b) combining said received signal samples with said channel tap estimates and 

4 said estimate of impairment correlation properties to produce Z(n) values; 

5 (b) combining said channel tap estimates with one another and said estimate of 

6 impairment correlation properties to produce S(k,n) values; and 

7 (c) combining said Z(n) and S(k,n) values with said hypothetical symbol 

8 sequences to produce branch metrics. 
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1 11. An apparatus for reducing the effects of signal fading, time dispersion, and 

2 Interference in a radio communications system, comprising: 

3 a digital transmitter for generating and transmitting a signal representing a 

4 transmitted symbol sequence; 

5 a receiver for receiving said signal on at least two separate antennas; 

6 an analog to digital converter for converting said signal into received signal 

7 samples for each of said antennas; 

8 a channel tap estimator for estimating channel taps for each of said antennas to 

9 produce channel tap estimates; 

10 an impairment correlation estimator for estimating impairment correlation 

1 1 properties to produce an estimate of impairment correlation properties; 

12 a branch metric processor for forming branch metrics using said received signal 

13 samples, said channel tap estimates, and said estimate of impairment correlation 

14 properties; and 

15 a sequence estimation processor employing said branch metrics to estimate 

16 said transmitted symbol sequence. 

1 12. The apparatus of claim 11, in which said analog to digital converter employs log- 

2 polar signal processing. 

1 13. The apparatus of claim 1 1 ( in which said branch metric processor comprises: 

2 a symbol sequence generator for generating hypothetical symbol sequences; 

3 means for filtering said hypothetical signal sequences with said channel tap 

4 estimates to produce hypothesized received signal samples for each antenna; 

5 means for subtracting said hypothesized received signal samples from the 

6 received signal samples to produce hypothesized error signals; and 

7 means for processing said hypothesized error signals with said estimate of 

8 impairment correlation properties to produce branch metrics. 
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1 14. The apparatus of claim 11 ( in which said impairment correlation estimator 

2 comprises: 

3 a symbol sequence generator for generating tentative detected symbol 

4 sequences; 

5 a digital filter for filtering said tentative detected symbol sequences with said 

6 channel tap estimates to produce detected signal samples for each antenna; 

7 means for subtracting said detected signal samples from the received signal 

8 samples to produce detected error signals; and 

9 means for processing said detected error signals with said estimate of 

10 impairment correlation properties to produce an update of the estimate the impairment. 

11 correlation properties. 

1 15. The apparatus of claim 1 1 , in which said branch metric processor comprises: 

2 a symbol sequence generator for generating hypothetical symbol sequences; 

3 means for calculating pre-computed values for all of said hypothetical symbol 

4 sequences using said channel tap estimates and said estimate of impairment 

5 correlation properties; and 

6 means for processing said received signal samples with said pre-computed 

7 values to produce branch metrics. 

1 16. The apparatus of claim 1 1 1 in which said branch metric processor comprises: 

2 a symbol sequence generator for generating hypothetical symbol sequences; 

3 means for combining the received signal samples with the channel tap 

4 estimates and said estimate of impairment correlation properties to produce metric 

5 multipliers e(j,n) and for combining said channel tap estimates with one another and 

6 said estimate of impairment correlation properties to produce metric multipliers f(j,n) 

7 and g(j,k,n); and for combining said metric multipliers with said hypothetical symbol 

8 sequences to produce branch metrics. 
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1 1 7. The apparatus of claim 1 1 1 in which said branch metric processor comprises: 

2 a symbol sequence generator for generating hypothetical symbol sequences; 

3 means for combining said received signal samples with said channel tap 

4 estimates and said estimate of impairment correlation properties to produce Z(n) 

5 values and combining said channel tap estimates with one another and said estimate of 
8 impairment correlation properties to produce S(k,n) values and for combining said Z(n) 
7 and S(k,n) values with said hypothetical symbol sequences to produce branch metrics. 

1 18. The apparatus of claim 11, in which said sequence estimation processor employs 

2 the Viterbi algorithm. 



1 19. A method for reducing the effects of signal fading, time dispersion, and interference 

2 in a radio communications system, comprising the steps of: 

3 (a) generating and transmitting a signal representing a transmitted symbol 

4 sequence; 

5 (b) receiving said signal on one antenna; 

6 (c) processing said signal to produce received signal samples for said antenna 

7 such that multiple signal samples per symbol period are produced; 

8 (d) estimating channel taps for each of said antennas to produce channel tap 

9 estimates; 

10 (e) estimating impairment correlation properties to produce an estimate of 

1 1 impairment correlation properties; 

12 (f) forming branch metrics in a branch metric processor using said received 

13 signal samples, said channel tap estimates, and said estimate of impairment correlation 

14 properties; and 

15 (g) employing said branch metrics in a sequence estimation algorithm to 

16 estimate said transmitted symbol sequence. 
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1 20. A method according to claim 19 ( in which said forming step (0 comprises: 

2 a) generating hypothetical symbol sequences; 

3 b) combining the received signal samples with channel tap estimates and said 

4 estimate of impairment correlation properties to produce metric multipliers e(j ( n); 

5 c) combining said channel tap estimates with one another and said estimate of 

6 impairment correlation properties to produce metric multipliers f(j,n) and g(j,k ( n); and 

7 d) combining said metric premultipliers with said hypothetical symbol sequences 
a to produce branch metrics. 



1 21 . A method according to claim 19, in which said forming step (0 comprises: 

2 a) generating hypothetical symbol sequences; 

3 b) combining the received signal samples with channel tap estimates and said 

4 estimate of impairment correlation properties to produce metric multipliers e(j,n); 

5 c) combining said channel tap estimates with one another and said estimate of 

6 impairment correlation properties to produce metric multipliers g(j,M); and 

7 d) combining said metric premultipliers with said hypothetical symbol sequences 

8 to produce branch metrics. 



1 22. A method according to claim 19, in which said forming step (0 comprises: 

2 a) generating hypothetical symbol sequences; 

3 b) combining said received signal samples with said channel tap estimates and 

4 said estimate of impairment correlation properties to produce Z(n) values; 

5 (b) combining said channel tap estimates with one another and said estimate of 

6 impairment correlation properties to produce S{k t n) values; and 

7 (c) combining said Z(n) and S(k,n) values with said hypothetical symbol 

8 sequences to produce branch metrics. 
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1 23. An apparatus for reducing the effects of signal fading, time dispersion, and 

2 Interference in a radio communications system, comprising: 

3 a digital transmitter for generating and transmitting a signal representing a 

4 transmitted symbol sequence; 

5 a receiver for receiving said signal on one antenna; 

6 an analog to digital converter for converting said signal into received signal 

7 samples for said antenna such that multiple signal samples per symbol are produced; 

6 a channel tap estimator for estimating channel taps for each of said antennas to 

9 produce channel tap estimates; 

10 an impairment correlation estimator for estimating impairment correlation 

11 properties to produce an estimate of impairment correlation properties; 

12 a branch metric processor for forming branch metrics using said received signal 

13 samples, said channel tap estimates, and said estimate of impairment correlation 

14 properties; and 

15 a sequence estimation processor employing said branch metrics to estimate 

16 said transmitted symbol sequence. 

1 24. The apparatus of claim 23, in which said branch metric processor comprises: 

2 a symbol sequence generator for generating hypothetical symbol sequences; 

3 means for combining the received signal samples with the channel tap 

4 estimates and said estimate of impairment correlation properties to produce metric 

5 multipliers e(j,n) and for combining said channel tap estimates with one another and 

6 said estimate of impairment correlation properties to produce metric multipliers f(j.n) 

7 and gQ,k,n); and for combining said metric multipliers with said hypothetical symbol 

8 sequences to produce branch metrics. 
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1 25. The apparatus of claim 23, in which said branch metric processor comprises: 

2 a symbol sequence generator for generating hypothetical symbol sequences; 

3 means for combining said received signal samples with said channel tap 

4 estimates and said estimate of impairment correlation properties to produce Z(n) 

5 values and combining said channel tap estimates with one another and said estimate of 

6 impairment correlation properties to produce S(k,n) values and for combining said Z(n) 

7 and S(k,n) values with said hypothetical symbol sequences to produce branch metrics. 
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